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1.  INTOODUCTON 


Metal  matrix  composites  (MMCs)  have  been  sought  as  a  strong  candidate  for  high 
performance  structural  materials  for  use  in  both  room  and  high  temperature 
environments.  The  thermo-mechanical  behavior  of  various  types  of  MMCs  has  recently 
been  reviewed.^^^  Although  the  mechanical  behavior  of  MMCs  in  room  temperature  use 
environments  was  superior  to  that  of  the  unreinforced  metal,  it  is  the  high  temperature 
use  environment  in  which  degradation  in  the  mechanical  behavior  of  the  MMCs  has  been 
recognized.  Typical  high  temperature  use  environments  that  MMCs  experience  are  creep 
and  thermal  cycling.  Normally  the  creep  rate  of  a  MMC  is  much  less  than  that  of  the 
unreinforced  metal,  thus  exhibiting  the  merit  of  reinforcing  the  metal  by  fibers.  In  the 
long  term,  although  the  interface  of  the  matrix  and  fiber  may  become  degraded  due  to 
diffusion  of  matrix  metal  and  fiber  materials.  It  is,  however,  thermal  cycling  that  often 
limits  the  use  of  otherwise  promising  MMCs  in  high  temperature  applications. 
Degradation  of  the  mechanical  properties  and  dimensional  instability  of  MMCs  subjected 
to  thermal  cycling  have  been  reported. Table  1  shows  a  summary  of  the  degradation 
modes  of  various  types  of  MMCs  reported  so  far  in  the  literature.^^^^  The  degradation 
modes  stated  in  Table  1  stem  from  the  mismatch  in  the  coefficient  of  thermal  expansion 
(CTE)  of  the  matrix  metal  and  fiber  (normally  made  of  refractory  metal  or  ceramics).  The 
CTE  mismatch  strain  is  equal  to  (af-am)AT  where  a(  and  are  CTE  of  the  fiber  and 
matrix  metal,  and  AT  is  the  temperature  difference  undergone  by  the  MMC.  This  CTE 
mismatch  strain  induces  an  internal  stress  field  with  its  maximum  in  and  around  a  fiber, 
resulting  in  the  yielding  of  the  matrix  if  its  magnitude  is  large.  Under  thermal  cycling, 
this  internal  stress  field  and  cyclic  plastic  straining  of  the  matrix  metal  with  alternative 
sign  is  equivalent  to  low  cycle  (mechanical)  fatigue  described  by  Manson^^^V  This 
thermal  cycling  fatique  resuits  in  the  microdamage  of  the  matrix-fiber  interface:  voids  in 
W/Cu  composite^^\  cracks  in  B/6061  Al  composite(^^  and  also  macrodamage  such  as 
dimensional  changes,  which  are  unacceptable  since  MMC  components  undergoing  such 
dimensional  changes  may  not  function  properly  in  a  use  environment. 
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Table  1.  A  Summary  of  Previously  Reported  Thermal  Cycling  Damage  of  MMCs 


MMC  Systems 

Thermal  Cycling  Conditions 

Degradation  Modes 

References 

B/1100,  2024  and  6061  Al 

Tmln  “  RT 

Tmax  “  3^5,  365, 

N  -  6000 

surface  roughening. 

Internal  notohing  of  one 
matrix 

(2) 

B/6061  Al 

Tmin  “  ST 

Tmax  -  455"C 

N  •  1000,  2059,  and 

cmcklng  in  the  interfacial 
brittle  zorte 

(4) 

B/1100,  6061  Al 

Tmin  “  50“C 

Tmax  -  320,  420*0 

reaction  product. 

AIB2  caused  flawed  B  fibers 

(5) 

a-Al203  (FP)/A1-U 

Tmln  - 

Tmax  ■  500-C 

N  -  205,  1440,  2830 

voids  at  the 

reaction  product  LigO  • 

5AI203 

(8) 

a-Al203(FP)/ 

EZ33Mg 

Tmln  -  50*C 

Tmax  -  250,350*0 

N  -  up  to  3000 

voids  at  the  interface, 
matrix  cracMng,  matrix 
softening 

(11) 

Carbon  (M40)/ 

5056  Al 

Tmln  ”  nT 

Tmax  -  250-400*C 

N  -  10-1000 

dabonding  al  the  inlerface 

(13) 

SICw/2024  Al 

Tmln  “  RT 

Tmax  -  400,500*C 

N  .  100,1000 

dimensional  change  in  three- 
directions 

(12) 

W/Cu 

Tmln  -  200,400*0 

Tmax  -  800*C 

N  •  up  to  4950 

siding  at  the  interface 

(6) 

W/Cu 

Tmln  -  200*C 

Tmax  ■  800*C 

N  -  200 

voids  at  the  interface 

(7) 

W-1%Th02/ 

FsCrAlY 

Tmln  -  29*C 

Tmax  -  1204->C 

N  -  1000 

surface  roughening 

(3).  (15) 

Mo-T2M/FsCrAIY 

Tmln  -  29*C 

Tmax  •  1093,1204*0 

N  -  100,1000 

surface  cracMng, 
recrystaNizalion  of  TZM 

(3) 

W/sup«rafloys 

Tmln  •  RT 

Tmax  -  up  to  1100*0 

cracMng  at  the  interface 

(9) 

W-2wl%Th02/  stainless 
steel  Kovar,  Inconel  718, 
Hastellog  X 

Tmln  -  90*0 

Tmax  -  «25’C 

N  -  up  to  3600 

cracMng  and  dabonding  of  the 
interface 

(18) 

W/FeCrAlY 

Tmln  “  RT 

Tmax  -  350,1000, 

1090*0 

voids  in  the  matrix 

(14) 

W-1%  Th02/ 

FeCrAlY 

Tmln  “  RT 

Tmax  -  1100*0 

cracks  In  interface 

(1«) 
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In  this  report,  an  attempt  is  made  to  study  experimentally  and  analytically  the 
mechanisms  of  dimensional  changes  of  thermally-cycled  MMCs.  To  this  end,  we  have 
chosen  as  a  target  MMC  a  W-Th02  fiber/FeCrAlY  matrix  (W/FeCrAlY)  composite.  The 
experimental  study  of  the  W/FeCrAlY-  composite  will  be  discussed  in  Section  2. 
Analytical  modeling  of  dimensional  change  will  be  summarized  in  Section  3.  A 
comparison  will  then  be  made  between  the  present  experimental  results  and  the 
predictions  based  on  the  proposed  analytical  model  in  Section  4.  The  conclusion  is  given 
in  Section  5. 

2.  EXPERIMENT 

2. 1  MMC  Specimens 

In  this  study,  we  have  used  four  different  types  of  W/FeCrAlY  composites;  aligned 
continuous  fiber  MMC  (CF),  aligned  short  fiber  MMC  of  fiber  aspect  ratio  5  (ASF-5)  and  20 
(ASF-20),  and  a  hybrid  MMC  (AF).  In  the  hybrid  MMC,  chopped  AI2O3  fibers  are  combined 
with  the  FeCrAlY  metal  matrix  of  the  ASF-20  composite  with  the  intention  of  modifying 
the  otherwise  large  CTE  mismatch  strain  that  would  exist  between  the  W  fiber  and 
FeCrAlY  matrix.  These  MMC  types  are  summarized  in  Table  2. 


Table  2.  MMC  types  used 


MMC  Type 

Name 

Vol.  Frac.  of  W 
fiber  (%) 

Vol.  Frac.  of 
AI2O3  fibers(%) 

aligned  continuous 
fiber  MMC 

CF 

31.9 

0 

aligned  short  fiber 
MMC  (aspect  ratio-5) 

ASF-5 

19.4 

0 

aligned  short  fiber 
MMC  (aspect  ratio-20) 

ASF-20 

19.5 

0 

hybrid  MMC 

AF 

18.4 

10.0 

Processing  of  the  MMC  specimens  are  stated  in  subsection  2.1,  and  thermal  cycling 
testing  is  described  in  subsection  2.2.  Finally  the  results  of  the  thermal  cycling  tests  are 
presented  in  subsection  2.3. 
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2.2  Processing  of  MMC  Specimens 


The  processing  route  of  the  MMC  specimens  involves  several  steps:  processing  of  a 
matrix  mat  through  tape  casting,  processing  of  fiber  mats  (dirferent  for  the  case  of 

continuous  fibers,  and  short  fibers),  burn-off  of  polymer  binders,  consolidation  by  hot 

^  '  * 

isostatic  press  (HIP),  and  machining  into  specimens.  The  processing  route  is  illustrated 
schematically  in  Fig.  1.  The  matrix  metal  FeCrAlY  was  supplied  by  General  Electric  Corp. 
in  powder  form  with  a  nominal  diameter  of  15  pm.  The  W-Th02  fibers  were  supplied  by 
General  Electric  in  a  wire  form  of  nominal  diameter  of  250  pm,  while  AI2O3  fibers  (FP 
trade  name)  were  supplied  by  DuPont  and  their  nominal  diameter  is  19  pm. 

In  the  case  of  a  continuous  fiber  MMC,  fiber  mats  are  cut  from  wound  fiber  sheets 
prepared  by  filament  winding  where  a  polymer  binder  was  used.  A  similar  procedure  was 
utilized  for  producing  the  short  fiber  MMCs,  although  alignment  of  the  short  fiber  mats 
required  additional  care. 

The  polymer  burn-off  operation  is  necessary  to  remove  the  various  polymer 
materials  used  in  making  composite  preforms.  The  most  desirable  process  is  the 
diffusion  and  reaction  of  atmospheric  oxygen  with  polymer  materials  producing  CO2  and 
H2O  gaseous  reaction  products,  while  minimizing  the  oxidation  of  the  metal  powder. 
After  a  trial  and  error  process,  we  found  that  the  best  burn-off  schedule  consists  of  a 
rapid  temperature  climb  to  200*’C  followed  by  a  slow  two  hour  ramp  to  300°C.  The 
temperature  was  held  at  SOO^C  for  two  hours  followed  by  ramp  from  300*’C  to  350*’C  for 
two  hours.  Finally  the  preform  was  cooled-down  to  room  temperature.  The  entire  burn- 
off  procedure  required  approximately  nine  hours. 

The  HIP  process  used  required  the  temperature  and  pressure  to  be  increased  to 
1150°C  and  1500  MPa,  respectively  in  two  hours.  Both  conditions  were  then  held  for  two 
hours,  followed  by  a  rapid  release  in  pressure  and  slow  cool  in  temperature  over  night  so 
as  to  allow  creep  and  diffusional  relaxation  processes  to  relax  the  internal  stresses  in 
the  MMC. 
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MMC  specimens  are  so  hard  conventional  machining  methods  are  impossible  to  apply 
because  they  are  too  time  consuming.  We  thus  employed  abrasive  waterjet  cutting  method 
to  rough  machine  the  HIP  bulk  specimen  into  rod  shapes  followed  by  precision  machining 
by  tungsten  carbide  cutting  tools.  If  W-ThOs  fibers  become  exposed  to  the  atmosphere 
due  to  machining,  faster  oxidation  of  the  MMC  specimen  results.  To  avoid  this,  a  thermal 
spray  of  metal  FeCrALY  powder  was  applied  to  the  surface  of  the  MMC  specimens,  with  the 
thickness  of  coating  being  nominally  0.2  mm.  Further  details  of  the  MMC  processing  have 
been  given  in  our  previous  report.^^®^ 

2.3  Thermal  Cycling  Testing  and  Dimensional  Changes 

The  apparatus  for  thermal  cycling  testing  (hereafter  thermal  cycler)  is  shown  in 
Fig.  2  indicating  two  furnaces.  Each  furnace  is  a  fluidized  bed  so  as  to  ensure  rapid  heat 
transfer  to  a  MMC  specimen.  The  specimen  shuttles  within  10  seconds  between  the  low 
temperature  (Tn^jn)  and  high  temperature  (T^ax)  furnaces  with  programmed  time  sequence 
controlled  by  a  DC  motor  time  switch.  Tn,ax  can  be  as  high  as  1400°C  while  Tn^in  can  be 
reduced  to  room  temperature.  Examples  of  a  typical  temperature-time  relation  of  a  MMC 
specimen  are  shown  in  Fig.  3  where  two  cases  are  shown;  with  air  cooling  (dashed  line) 
and  without  air  cooling  (solid  line).  In  the  present  study,  Tmax  'S  set  equal  to  1100°C, 
while  three  different  T„,jn  are  used:  352,  534,  and  SOOX. 

The  longitudinal  (axial)  length  of  the  MMC  specimen  is  measured  before  and  after 
thermal  cycles  with  a  needlepoint  micrometer.  The  resolution  of  the  micrometer  is 
approximately  20  pm.  With  a  specimen  gauge  length  of  28  mm,  the  smallest  resolvable 
difference  in  the  data  is  approximately  0.001.  A  summary  of  longitudinal  (axial) 
dimensional  change  (A 7”)  where  a”  and  "  are  the  change  in  length,  and  initial  length  of  a 
specimen,  respectively  is  given  in  Table  3.  The  test  specimen  ID  number  in  Table  3  is 
defined  such  that  the  first  letters  stands  for  type  of  MMC,  the  first  number  for  T^ax*  aid 
the  second  number  for  T^in-  Some  of  MMC  specimens  tested  under  severe  conditions 
exhibited  excess  damage,  thus  those  specimens  are  not  usable  for  measurement  of 
dimensional  changes  (marked  by  *  in  Table  3). 
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Rgure  2.  Schematic  Diagram  of  the  Metallurgical  Thermal  Cycling  Machine. 


Table  3.  Longitudinal  Dimensional  Change  (Strain) 


Test  Specimen  ID 
Number 

Continuous 

CF-1 100-352 

.004683 

.02442 

• 

fiber  MMC 

CF-1 100-534 

.002257 

.00755 

* 

CF-1 100-800 

— 

damaged 

* 

ASF-20-1100-352 

.004252 

.01494 

* 

Aligned  short 

ASF-20-1100-352 

— 

.01577 

fiber  MMC 

ASF-20-1100-534 

.01080 

.02512 

* 

ASF-20-1100-800 

— 

.007162 

* 

ASF-5-1100-352 

-.000644 

.000270 

.00670 

ASF-5-1100-534 

.000911 

.002523 

.01514 

ASF-5-1100-800 

.003892 

— 

Hybrid  composite 

AF-1 100-352 

-.001606 

-0.000966 

.00734 

AF-1 100-534 

-.000588 

.002151 

.01203 

AF-1 100-800 

.000535 

— 

Those  specimens  marked  by  -  sign  were  not  tested  due  to  limited  number  of  specimens. 
The  results  of  the  axial  dimensional  change  of  the  four  types  of  MMC  specimens  with  Tmin 
=  352*0  and  T^ax  =  1100*C  are  plotted  as  a  function  of  the  number  of  thermal  cycles  (N)  in 
Fig.  4.  Figure  4  indicates  that  the  increasing  order  of  axial  dimensional  changes  observed 
is  continuous  fiber  MMCs  (CF),  aligned  short  fiber  MMCs  with  fiber  aspect  ratio  20  (ASF- 
20),  and  with  5  (ASF-5),  and  hybrid  composite  which  is  a  combination  of  ASF-20  and  10% 
volume  chopped  AI2O3  fibers  in  the  matrix  (AF).  This  order  is  expected  since  the  CTE 
mismatch  strain  is  larger  for  CF  than  for  both  ASF-5  and  AF.  Table  2  also  shows  their 
volume  fractions,  CF  having  31.9%  (largest)  and  the  other  having  almost  the  same  volume 
fraction  (»  19%).  A  comparison  between  the  experimental  and  analytical  results  of 
dimensional  change  will  be  given  in  Section  4. 


2.4  Observation  of  Damages 


In  order  to  assess  the  extent  of  the  damages  induced  in  the  thermally  cycled  MMCs, 
we  observed  both  uncycled  and  cycled  MMC  specimens  by  optical  microscope,  scanning 
electron  microscope  (SEM)  and  transmission  electron  microscope  (TEM). 
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(a)  (b) 


Figure  5.  Scanning  Electron  Micrograptis  of  Hybrid  Composite  Microstructure,  (a)  As 
HIPed,  (b)  After  500,  1100-352  Thermal  Cycles. 
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Figure  5  shows  a  normal  view  of  a  hybrid  MMC  (AF)  specimen  obtained  by  an  optical 
microscope  where  (a)  and  (b),  denote  the  as-received  and  500  cycled  specimens, 
respectively.  It  is  noted  in  Fig.  5  that  the  dark  islands  of  elliptical  shape  are  the  chopped 
AI2O3  fibers  mixed  with  the  matrix.  All  the  thermal  cycled  MMC  specimens  were 
subjected  to  the  same  thermal  cycling  conditions;  T^ax  =  HOOX,  Tmin  =  352°C.  and  with 
high  temperature  duration  of  x  =  4  minutes.  From  the  micrographs  taken  from  the  as- 
processed  and  thermal  cycled  MMCs,  one  can  characterize  the  micro-damage  associated 
with  the  as-received  and  as-thermally  cycled  MMCs.  The  thickness  of  circumferential 
reaction  zone  and  crack  density  (number  of  cracks  per  fiber)  within  the  zone  are  measured 
as  a  function  of  number  of  thermal  cycles,  and  the  results  are  plotted  in  Fig.  6(a)  and  (b). 
respectively.  Figure  6  indicates  that  hybrid  composite  (AF)  exhibited  the  best 
performance  among  the  others  whereas  CF  and  ASF-20  give  rise  to  larger  crack  density 
and  thicker  reaction  zone. 

A  limited  amount  of  TEM  study  was  conducted  on  the  thermal  cycled  ASF-20  and  AF 
specimens^^^)  with  the  aim  of  determining  the  effect  of  AI2O3  fibers  added  to  AF 
specimen  in  terms  of  dislocation  morphology.  Figures  7(a)  and  (b)  show  TEM  photographs 
of  the  dislocation  morphology  of  the  matrix-reaction  zone  interface  in  the  ASF-20 
specimen  and  that  of  the  matrix-Ai203  fiber  interface  in  the  AF  specimens,  respectively. 
These  TEM  photos  suggest  that  the  addition  of  AI2O3  chopped  fibers  into  the  matrix  of  the 
AF  (hybrid)  specimen  relaxes  the  otherwise  high  stress  concentration  at  the  matrix-W- 
Th02  fiber  interface  induced  by  the  punching  of  dislocations  [20]  by  the  CTE  mismatch 
strain.  At  the  same  time,  the  matrix  with  AI2O3  fibers  becomes  strengthened.  These 
mechanisms  contribute  to  the  increase  in  thermal  cycling  resistance  of  hybrid  MMC. 

In  order  to  obtain  some  clues  to  the  accelerated  rate  of  dimensional  changes 
observed  in  the  thermally  cycled  MMCs  with  higher  number  of  cycles.  Fig.  8,  we  conducted 
SEM  observations  on  longitudinal  sections  of  the  cycled  specimens,  focusing  on  the 
matrix-fiber  (W-Th02)  interface.  To  this  end.  numerous  SEM  photographs  were  taken 
along  the  interfacial  zones  of  ASF-5,  ASF-20,  and  AF  specimens.  Figures  9(a)  and  (b) 
show  the  SEM  micrographs  of  the  interface  near  the  fiber-end  of  the  as-processed  and 
thermal-cycled  N-500)  ASF-5  specimen,  respectively.  It  Is  clearly  observed  in  Fig.  9  that 
thermal  cycling  to  N  -  500  resulted  in  debonding  of  the  matrix-fiber  interface.  In  order  to 
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Figure  7.  Transmission  Electron  Micrographs  of  Thermal  Cycled  Composite,  (a)  W(Fe,CR) 
Reaction  Zone  -FeCrAlY  Matrix  Interface,  (b)  Al203-FeCrAIY  Interface. 
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Dimensional  Change  as  a  Function  of  Thermal  Cycles,  Non>Unear  for  the 
Second  Stage. 
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Rgure  9.  SEM  Photographs  Showing  Thermal  Cycling  Induced  Interfacial  Damage  and 
Results  of  Damage  Quantification  by  Image  Analysis. 


quantify  the  damage  near  the  interface,  the  SEM  photographs  enhanced  to  provide  black  and 
white  pictures  where  black  domains  denote  voided  (or  damaged)  area.  The  enhanced 
pictures  were  scanned  by  an  Apple  scanner  connected  to  Macintosh-ll  computer  system  to 
create  a  scan  file  where  each  black  pixel  represented  0.0675  square  micron  of  damage.  By 
using  computer  software  developed  to  count  the  number  of  black  pixels  in  the  scanfile,  we 
converted  the  damage  profile  shown  in  the  SEM  photographs  to  an  area  distribution  of 
damage  along  the  fiber  length.  Examples  of  the  area  distribution  of  damage  are  shown  in 
Fig.  9(c)  where  the  dashed  and  solid  lines  denote  the  damage  area  of  the  500  and  750 
cycled  ASF-5  specimens,  respectively.  Figure  10  shows  a  comparison  of  the  damage  area 
plots  of  ASF-20  and  hybrid  MMCs,  both  thermal  cycled  to  N  «  500.  Figures  9  and  10 

indicate  that  for  N  =  500,  the  hybrid  composite  gives  rise  to  the  smallest  amount 

interfacial  damage  while  the  ASF-20  specimen  yields  the  largest  amount  of  damage.  Thus 
for  the  same  fiber  aspect  ratio  {=  20),  the  thermal  cycling  resistance  of  the  hybrid 
composite  improved  over  that  of  the  non-hybrid  composite  (ASF-20).  It  follows  from 
Figs.  9  and  10  that  the  extent  of  the  damaged  area  is  the  largest  for  the  ASF-20 

specimens  and  is  the  smallest  for  the  AF  specimen.  The  thermal  cycling  damage  appears 

to  start  from  fiber-ends  and  progressively  extends  towards  the  center  of  fiber  as  the 
number  of  thermal  cycles  increases.  This  is  illustrated  schematically  in  Fig.  11.  This 
interfacial  damage  reduces  the  load  transfer  capacity  from  the  matrix  to  the  short  fiber 
which  is  a  key  strengthening  mechanism  for  short  fiber  MMCs.  The  damage  would  also 
have  a  strong  effect  on  the  dimensional  change. 

3.  ANALYTICAL  MODEL 

The  present  analytical  model  is  aimed  at  predicting  dimensional  changes  induced  in 
thermally  cycled  MMCs.  Typical  experimental  data  of  axial  dimensional  change  (el)  of 
several  MMCs  is  summarized  by  Fig.  8:  The  el  -  N  curve  consists  of  basically  two  stages, 
the  first  stage  is  where  the  el  -  N  relation  is  linear,  the  second  stage  where  it  becomes 
highly  non-linear  indicating  progressive  damage.  The  first  stage  corresponds  to  the  case 
of  perfectly  bonded  interfaces,  while  the  second  stage  to  the  case  of  progressively 
debonded  interfaces.  We  shall  treat  these  two  cases  separately  in  the  following.  Most 
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Figure  10.  Interfacial  Damage  Quantification  by  Image  Analysis  Comparing 
Cycled  Non-Hybrid  and  Hybrid.  Aspect  Ratio  20  Composite  Materials. 
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of  the  formulation  used  in  the  analytical  model  for  the  first  stage,  however,  can  also  be 
used  for  the  second-stage  model. 

3.1  Analytical  Model  for  Perfectly  Bonded  Interface 

The  temperature-time  history  for  thermal  cycling  is  idealized  as  shown  in  Fig.  12 
where  the  specimen  is  initially  subjected  to  Tmax  at  which  it  is  assumed  stress-free,  and 
then  experiences  the  temperature  sequence:  a-»b-4c->d->e. 

The  main  assumptions  used  in  the  present  model  are: 

(1)  For  the  cooling  process  (a  -¥  b),  Tmax  -  Tmin  is  large  enough  to  yield  the  matrix 
metal  of  the  MMC  with  yield  stress  (oy)  at  some  temperature  (Tyieid)  while  the 
fiber  remains  elastic  throughout  the  cycle. 

HT 

(2)  At  Tmax.  the  matrix  can  yield  with  the  high  temperature  yield  stress  (o  y  )  and 
also  can  creep  according  to  a  power-law  of  the  Oorn  type.f^^) 

(3)  For  the  majority  of  the  analysis,  short  fibers  are  aligned  in  one  direction.  Fig. 
13.  We  will  discuss  later  the  effect  of  fiber  misorientation  on  dimensional 
change.  Fig.  14. 

(4)  Both  the  matrix  and  fiber  are  isotropic  in  stiffness  (E  .  v)  and  CTE(a). 

Cooling  Process  a  b 

Referred  to  Fig.  12,  as  the  temperature  drops  from  Tmax  to  Tmin.  the  matrix  metal 

undergoes  two  different  deformations:  elastic  (from  Tmax  fyiaid)  and  plastic  (Tyi^id 

fmin)*  Hence,  we  will  first  determine  Ty|,id  at  which  the  matrix  starts  yielding.  The  CTE 
mismatch  strain  that  exists  in  the  domain  of  fiber,  a*  can  be  given  by 


a*  -  (Of  -  am,  ocf  -  Om,  Sf  -  Om  .  0,  0,  0)  AT 


(1) 


T 

max 


e 


T 

yield,  heating 
T 

yield,  cooling 

T  . 
mm 


a  d 

- ► 

b  c 


time 


Figure  12.  Idealized  Time  (t)  versus  Temperature  (T)  Relation  Used  in  the  Present  Model. 


Rgure  13.  Domain  of  the  Analytical  Model 
for  an  Aligned  Short  Fiber  MMC 
and  Definitions  of  Coordinate 
System  Used. 


Figure  14.  Domain  of  the  Analytical 

Model  for  a  Misoriented  Short 
Fiber  Composite  and  Definition 
of  Coordinate  System  Used. 
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where  the  components  of  a*  are  11,  22.  33.  23.  31,  and  12,  from  the  right  to  left  in  order. 

In  eq.  (1)  af  and  am  are  the  CTEs  of  the  fiber  and  matrix,  respectively,  and  AT  is  the' 
temperature  change  (equal  to  Tyjei,^  -  Tmax>  i'l  case  of  the  cooling  process  a  y).  The 
volume  averaged  stress  in  the  fiber  and  matrix.  Om  are  calculated  using  the  Eshelby’s 

equivalent  inclusion  method,  (30)  see  Appendix  A,  where  the  eigenstrain  ej  is  replaced  by 
a*,  Eq.(1). 

oj  -  (1-f)  Cm  •  (S  •  e*  -  e* )  (2) 

Om  =  Cm  (S  •  f  *  -  e* )  (3) 


where  Cm  is  the  stiffness  tensor  of  the  matrix  (Young’s  modulus  Em  and  Poisson's  ratio 
vm  for  isotropic  matrix),  S  is  Eshelby’s  tensor  which  is  a  function  of  the  fiber  aspect 
ratio  (7d)  and  vm  whose  explicit  expression  has  been  given  elsewhere(^).  In  eqs.  (1)  and 
(2)  e*  is  the  fictitious  eigenstrain  introduced  into  the  fiber  domain  (see  Appendix  A).  The 

A* 

dots  used  in  Eqs.  (2)  and  (3)  are  inner  product  between  tensors,  and  used  throughout  this 
paper.  For  the  matrix  to  yield  at  Tyj^id,  Om  must  satisfy  the  yield  condition: 


3  y’  y'  2 
2 


(4) 


where  Oy  is  the  matrix  yield  stress  at  Tyi^id  (set  equal  to  Tmin  for  simplicity  in  the 
analysis),  and  am^'  is  the  deviatoric  component  of  Omy*-  Eq.  (4)  provides  the  value  of  Tyi^id. 
As  the  temperature  is  further  reduced  from  Tyj,id  to  Tmfn.  the  corresponding  CTE  mismatch 
strain  a*  will  be  used  to  generate  plastic  strain,  ef*  at  Tmin  with  its  components  given  by 

»  (-«p/2.  -«p/2,  ep  0,  0,  0)  (5) 

where  ep  denotes  the  axial  component  of  the  plastic  strain  (i.e.,  along  the  xa-axis.  Fig. 
13),  which  is  still  unknown  at  this  stage  of  the  formulation,  but  will  be  determined  as 
follows.  A  given  plastic  strain  in  the  matrix  induces  the  internal  stress  in  the  matrix 
(Om)  and  in  the  fiber  (of)  which  can  be  calculated  by  the  Eshelby’s  equivalent  inclusion 
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method,  as  shown  in  Appendix  A  where  zj  •  -ep-  The  average  stress  in  the  matrix  at  Tf„i„ 
is  the  sum  of  and  o^,  the  latter  being  the  thermo-elastic  stress  obtained  for  ey  * 
with  AT  »  Tmin  -Tn,ax-  Assuming  a  perfectly  plastic  matrix,  this  average  matrix  stress  at 
Tniin  (its  axial  and  transverse  components  being  0^33  and  on.  respectively)  must  satisfy 
the  yield  condition: 


Om33  ~  ®m1 1  “  (6) 

Eq.  (6)  then  provides  the  value  of  ep,  which  will  contribute  to  the  macroscopic  plastic 
strain  of  the  MMC  at  T^in  (stage  b). 

ef  =  (1-f)ep 

(7) 

ef  -  -(1-f)ep/2 

where  the  subscripts  L  and  T  denote  the  component  along  the  longitudinal  (axial,  or  X3) 
and  transverse,  or  x^-,  X2-)  axes,  respectively. 

Heating  process  c  d 

As  the  temperature  changes  rapidly  from  T^in  (stage  c)  to  Tmax  (stage  d),  the  CTE 

mismatch  strain  defined  by  Eq.  (1)  with  AT  -  T„,ax  -Tmin  induces  the  internal  stress  in  the 

p  T 

fiber  and  matrix.  Hence  the  average  stress  in  the  matrix  at  stage  d  becomes  ,  since  o^ 

A* 

cancels  with  that  induced  during  the  cooling  process  (c  d).  The  matrix  average  stress 

p 

0^  must  satisfy  the  yield  condition  at  T^ax  effective  stress  exceeds  the  high 

H  T 

temperature  yield  stress  a  y  .  This  can  be  accomplished  by  introducing  an  additional 
plastic  strain  in  the  matrix, 

UTD  HTP 

The  average  matrix  stress  induced  by  (denoted  by  )  can  be  calculated  by  Eq. 

(A1)  where  ej  is  replaced  by  Then,  the  total  averaged  stress  in  the  matrix  at  stage 

p  fTrp  ~ 

d,  is  which  must  satisfy  the  yield  condition: 
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(8) 


,  P  HTP,  ,  P  HTP,  HT 
V®m33  +®m33)  “  l®mii  +®mil)  =~Oy 


It  should  be  noted  here  that  the  effective  (flow)  matrix  stress  at  stage  d  is  compressive. 
Additional  plastic  strain  at  stage  d,  contributes  to  the  macroscopic  plastic  strain, 

UTP  HTP  HTP 

(longitudinal  composite  denoted  by  and  transverse  component  by  e-]-  ), 
e”’’’’’  =  (1-f)eHTP 

(9) 

=  -(l-f)eHTp/2 

UT 

At  the  end  of  stage  d,  the  matrix  flow  stress  is  equal  to  o  y  . 


High  temperature  plateau:  d  ->  e 


During  the  isothermal  exposure  at  Tmax>  matrix  metal  stress  creeps  to  reduce  its 
flow  stress.  In  order  to  simulate  the  creep  relaxation  process  of  the  matrix,  we  shall 
subdivide  the  duration  of  the  high  temperature  plateau  (x)  into  m  time  increments.  For 
each  time  increment  At,  the  following  Dorn  type  creep  law  of  Dorn  type^^^^  is  applied  to 
calculate  the  resulting  creep  strain  increment. 


Aec-  - 


ronif  GbD, 

f  Qy 

rUJ 

^maxJJ 

At 


(10) 


where 

Aec  -  incremental  axial  creep  strain 
A  a  Dorn  constant 

Ont  -  matrix  flow  stress 

G  -  matrix  shear  modulus 

b  a  matrix  Burgers  vector 
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k  >  Boltzman's  constant 

0^  >  pre-exponential  constant  for  self-diffusion 

R  «  gas  constant 

n  -  power  law  constant 

Qy  =<  activation  energy 

The  minus  sign  on  the  right  hand  side  of  Eq.  (10)  indicates  that  the  matrix  is  in 
compression  during  the  high  temperature  plateau.  Once  the  creep  strain  increment  Aec  is 
obtained  from  Eq.  (10),  we  will  calculate  the  stress  induced  by  Aec,  by  Eq.  (A2) 
where 

ej  =  Ae® 


and  where 

Ae®  -  (-Aec/2,  -Aec/2,  Aec,  0,  0,  0)  (11) 

Then,  the  total  averaged  stress  in  the  matrix  becomes  Om  +  Aom,  from  which  the  matrix 
flow  stress  in  the  next  incremental  step.  is  calculated  as 

/ 

Om  “  (®m33  +  A®m33)  ”  (®m1 1  +^®m1l)  (^2) 

will  then  be  substituted  into  Eq.  (10)  to  calculate  the  next  incremental  creep  strain. 

At  the  end  of  the  incremental  process,  i.e.,  at  stage  e,  the  total  creep  strain  along  the  X3- 
axis  becomes 

ec-lASc  (13) 

ui 

Sc  contributes  to  the  macroscopic  plastic  strain  of  the  MMC,  e®  with  its  longitudinal 
and  transverse  component  ej  given  by 
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(1-f)ec 


c 

= 

(14) 

4  =  -(1-f)ec/2 

Dimensional  Change  per  Cycle.  Ae 

At  the  completion  of  one  thermal  cycle,  i.e.,  at  stage  e,  the  total  plastic  strain  Ae 
(dimensional  change  per  cycle)  can  be  calculated  from  Eqs.  (7),  (9),  and  (14). 

Along  the  fiber  axis; 

p  htp  c  , 

A£l  =  +  (15a) 


Along  the  transverse  direction: 
,  p  HTP  c 


(15b) 


Eq.  (15)  provides  the  dimensional  change  of  an  MMC  for  the  first  stage  of  el  -  N  curve 
where  the  matrix-fiber  interface  is  presumably  well  bonded. 

3.2  Analytical  Model  for  Progressively  Debonded  Interface 

The  formulation  derived  for  the  case  of  perfectly  bonded  interfaces  can  be 
applicable  to  that  of  progressive  debonded  interface  with  minor  modification.  The 
modification  is  to  account  for  the  effective  volume  fraction  of  fiber  fg  for  the  composite 
with  debonded  interfaces,  fe  depends  on  the  extent  of  the  debonded  interfaces.  Since  the 
stress  concentration  is  highest  at  fiber-ends,  the  debonding  is  more  likely  to  take  place 
there  first,  and  then  extends  toward  the  fiber  center  as  evidenced  by  SEM  photographs, 
Fig.  9. 


During  the  cooling  process  (a  ->  b.  Fig.  12),  the  radial  stress  {a„)  induced  at  the 
matrix-fiber  interface  due  to  GTE  mismatch  is  compressive,  while  for  the  heating  period 
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(c  -*  d),  Orr  at  the  interface  becomes  compressive.^^®^  This  leads  to  the  length  of  the 
debonded  interface  (d”)  for  cooling  being  smaller  than  that  for  heating.  The  matrix 
material  adjacent  to  the  tip  of  the  debonding  is  subjected  to  the  highest  shear  stress  (or 
strain,  Ay)  which  repeats  between  opposite  signs.  Therefore,  the  location  of  the  tip  at 
stage  a  (the  beginning  of  each  cycle)  will  extend  progressively  with  accumulated  number 
of  cycles  N  due  to  the  propagation  law  of  low  cycle  fatigue  given  by^^^) 

ay  =  yf(2N)“  (16) 

where  yf  >s  the  fracture  shear  strain  in  a  static  shear  test  (2N  1)  and  a  is  material 

constant.  The  sequence  of  progressive  debonding  interfaces  is  schematically  shown  in 
Fig.  15  where  the  cycles  (N-1)  -4  N  (N+1)  are  considered.  The  lengths  of  the  debonded 
interface  during  the  cooling  (a  b)  and  heating  (c  d  e)  at  N  cycles  are  denoted  by  d”]^ 

and  d”|^,  respectively. 

The  superscripts  (-•-)  and  (-)  corresponds  to  the  radial  stress  at  the  interface  being 
in  tension  and  compression,  respectively.  It  is  noted  in  Fig.  15  that  d”^^^  >  d”^  >  d”^.^, 
indicating  that  the  debonding  progresses  with  N.  In  reality,  this  is  governed  by  Eq.  (16), 
thus  the  extension  may  take  place  only  after  a  critical  number  of  accumulated  cycles.  Due 
to  the  unavailability  of  the  low  cycle  data  (Eq.  (16)),  we  employ  here  the  assumption  that 
the  extension  is  a  continuous  process  with  N.  The  change  in  the  length  of  the  debonded 
interface  d”^  'S  directly  related  to  that  of  the  effective  fiber  length  ”n,  which  in  turn  can 
be  interpreted  as  that  in  the  effective  volume  fraction  of  fiber,  f|g. 

It  is  clear  from  Eqs.  (7),  (9),  (14),  and  (15)  that  the  accumulated  inelastic  strain  per 
cycle,  Ae  is  dependent  on  (1-f)  where  f  is  the  effective  volume  fraction  of  fiber,  thus 
equal  to  f^j  for  a  given  cycle  N.  It  should  be  noted  again  that  f^  is  not  constant  during  this 
cycle  N:  it  is  larger  for  the  cooling  period  than  for  the  heating  period  and  high  temperature 
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Figure  15.  Progressive  Debonding  of  Fiber  Interfaces  in  a  Thermal  Cycled  Composite. 


plateau,  Fig.  15.  With  this  in  mind,  one  can  express  the  longitudinal  dimensional  change 
per  cycle  during  the  Nth  cycle  as: 


=  (■'“fN)  ®p  +  (^~^n)  (®htp  +  ec) 


(17) 


where  and  are  the  effective  volume  fractions  during  the  Nth  cycle  during  the  cooling 
(a  b)  and  heating  and  high  temperature  plateau  (c  -4  d  e),  respectively. 

In  actual  computation,  one  must  specify  the  rate  of  change  of  effective  fiber  volume 
fraction  between  N  and  N+i  cycles  focusing  on  the  compression  (a-b)  period,  dfN  is  defined 
as 


dlN  -  1n+i  -  In 


(18) 


One  also  needs  to  specify  the  differences  in  f^  between  the  heating  (c-d-»  e)  and  cooling 
(a-»b)  periods,  AfN  given  by 


AfN  -1n-1n 


(19) 


With  the  values  of  df^  and  Af^  given,  one  can  calculate  the  dimensional  change  for  the  N- 
th  cycle,  Ae^.  Then  the  accumulated  dimensional  change  up  to  N  cycles,  is  given  by 

N 

e^-Xae'  (20) 

1-1 

4.  DISCUSSION 

To  ensure  the  validity  of  the  present  model,  we  shall  compare  predictions  by  the 
present  model  with  those  of  existing  models,  with  the  experimental  results  reported  by 
others,  and  with  the  present  experimental  results.  We  will  discuss  the  predictions  for 
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the  first  stage  of  e-N  curve  in  Subsection  4.1,  followed  by  those  for  the  second  stage  in 
Subsection  4.2. 

4. 1  Dimensional  Change  for  the  First  Stage 


The  other  models  with  which  we  shall  compare  results  from  the  present  one  are 
Garmong’s  model, the  Yoda  et  al.  model, and  Taya-Mori  model.^^®^  For  the  purpose 
of  comparison,  we  use  the  data  of  a  hypothetical  MMC  (continuous  fiber)  given  by  Garmong: 


GTE  mismatch; 
fiber  volume  fraction: 
matrix  Young’s  modulus: 
matrix  shear  modulus: 
fiber  Young’s  modulus: 
matrix  yield  stress: 
constants  of  matrix  flow 
stress-plastic  strain  relations*: 

creep  constants** 


Thermal  cycling  condition: 


ttm-af  -  10  X  10"®/K 
f  -  0.5 

Em  “  6.9  X  lO^MPa 
Gm  =  2.62x  lO^MPa 
Ef«2.1  X  10®MPa 
Oy  -  0.56  -  0.07  X  T  (MPa) 

K*  =  690MPa 
n  .  0.5 

DoAGb/k  =  3.1  X  10“/sec 
Q/R  =  17600K 
S  »  4.4 

Tmax  »  700K 
Tmin  =  300K 
X  a  60  sec 


where  *  and  **  denote 

tx*  n 

Ofti  “  K  ®p 


and  the  creep  law  of  Dorn  type  defined  by  Eq.  (10),  respectively. 

With  regard  to  Garmong’s  model,  we  shall  use  the  model  corrected  by  Tyson.(^7) 
Analytical  results  of  dimensional  change  by  the  Yoda  et  al.  model  cannot  be  obtained  by 
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straight  application  of  the  Garmong  data,  thus,  we  use  here  the  upperbound  of  the 
dimensional  change  by  the  Yoda  et  al.  model.^®) 


HT 


1  -f 
f 


(21) 


Where  Oy  is  the  matrix  yield  stress  at  results  of  dimensional  change  per  cycle, 

Aei,.  of  the  hypothetical  MMC  are  plotted  as  a  function  of  fiber  volume  fraction,  f,  in  Fig. 

16  where  (a)  and  (b)  cover  the  entire  range  of  f,  and  the  range  of  f  focused  on  the  Ac|.  peak 

values  (f  »  0  ~  0.1),  respectively.  In  the  figure,  the  results  of  the  Garmong  model 

multiplied  by  (l-f)  are  also  shown,  which  becomes  zero  as  f  approaches  1.  It  is  noted  in 

UB 

Fig.  1 6  that  the  results  of  Ae  by  Yoda  et  al.  model  lie  far  above  those  of  the  other  models 

and  do  not  possess  a  peak  like  the  others.  Moreover  it  does  not  converge  to  zero  as  f  1 . 
It  also  follows  from  Fig.  16  that  the  dimensional  change  per  cycle  by  both  Garmong  and 
Taya-Mori  models  do  not  approach  to  zero  as  f  1.  The  present  model  exhibits  the 
correct  tendency;  Ael  ^  0  as  f  ^  1 ,  indicating  that  the  present  model  can  account  for  the 
interaction  between  fibers  at  higher  volume  fractions.  An  important  implication  of  Fig. 
16  is  that  there  exists  a  peak  value  of  ael  corresponding  to  a  small  volume  fraction  of 
fiber,  f  =  0.06  ~  0.065.  It  would  be  wise  then  for  us  to  avoid  this  critical  volume  fraction 
to  minimize  the  dimensional  change.  The  effect  of  fiber  aspect  ratio  predicted  by  the 
present  model  for  the  Garmong's  data  with  f  >  0.2  is  examined  and  the  results  are  shown 
in  Fig.  17,  suggesting  several  important  conclusions:  (1)  for  aspect  ratios  (7d)  larger 
than  6,  Ael  is  not  sensitive  to  7d,  (2)  Aei.  has  a  peak  at  about  7d  -  1.6,  and  (3)  a  particle 
MMC  where  the  reinforcement  aspect  ratio  is  1,  would  not  give  rise  to  any  dimensional 
change.  The  last  conclusion,  however,  does  not  lead  to  the  conclusion  that  particle  MMC  is 
free  of  thermal  cycling  damage;  instead,  the  cyclic  stress  in  an  around  a  particle  due  to 
thermal  cycling  may  result  in  the  debonding  of  the  matrix-particle  interface. 


A  limited  amount  of  experimental  data  are  available  in  the  literature  on  dimensional 
change  of  thermally  cycled  MMCs  during  the  first  stage  (Figure  8).  These  studies  include 
the  investigation  of  W/Cu  by  Yoda  et  al.(^^  SiC^/AI  by  Patterson  and  Taya.^^^^  and  the 
present  experimental  study  on  W-1%Th02/FeCrAiY.  The  input  data  of  these  MMC  systems 
for  the  present  model  are  summarized  in  Table  4. 
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Axial  Strain  Par  Tharmal  Cyela 


(a) 


(b) 


Figure  16.  Comparison  of  Existing  Thermal  Cycling  Deformation  Modeis  with  the  Present 
Model,  (a)  For  Entire  Range  of  Volume  Fraction,  (b)  Near  the  Peak  Region. 


Aapael  Ratio 


Figure  17.  Longitudinal  Strain  Per  Thermal  Cycle  (Acl)  Versus  Fiber  Aspect  Ratio  for  a 
Volume  Fraction  0.20  W-1%Th02  Reinforced  FeCrAlY  Composite. 
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Table  4.  Input-Data  for  Three  MMC  Systems 


W-Cu 

SiC-AI 

W-Th02-FeCrYAI 

Matrix 

En,(MPa) 

1 .1x10’® 

6.76x10* 

2.01x10® 

Vm 

0.35 

0.33 

0.33 

Om(10-®K-1) 

20.0 

24.7 

26.1 

Oy(MPa) 

68.9 

50.8 

220 

Oy^^''’(MPa) 

11.4 

34.5 

6.9 

Creep  Constant 

1.04x10^® 

3.41x10^2 

1.66x10^® 

s 

4.8 

5.08 

6.09 

Fiber 

Ef(MPa) 

3.58x10® 

4.27x10® 

3.58x10® 

Vf 

0.285 

0.17 

0.285 

af(10-®K'^) 

4.75 

4.3 

4.75 

Aspect  Ratio 

oo 

4.0 

5 

f 

0.1  and  0.5 

0.15 

0.2  and  0.3 

Thermal  Cycle 

^max(*’C) 

800 

400 

1100 

conditions 

Tmtn("C) 

200 

20 

352  and  534 

x(s) 

250 

125 

260 

The  prediction  based  on  the  present  model  for  the  axial  dimensional  change  per  cycle 
Aet  of  W/Cu  composite,  shown  by  a  solid  curve  in  Fig.  18,  is  compared  with  the 
experimental  results  of  Yoda  et  al.,  resulting  in  good  agreement.  In  the  case  of  SiCw/AI 
composite,  the  SiC  whiskers  are  distributed  in-plane  with  the  whisker  orientation  being 
semi-random,  i.e.,  not  aligned  in  a  specified  direction,  this  represents  a  departure  from 
the  present  model.  Fig.  13.  Hence,  the  dimensional  change  per  cycle  along  the  transverse 
direction,  Aej  is  focused  on  for  comparison,  for  the  dimensional  change  along  the 
thickness  direction  should  not  be  influenced  significantly  by  the  misorientation  of  fibers. 
Fig.  19  shows  the  prediction  of  Aej  by  the  present  model  (shown  by  solid  curve)  along  with 
the  experimental  result  of  Patterson  and  Taya  (filled  circle).  Again,  the  present  model 
agrees  with  the  experiment. 

In  order  to  apply  the  present  model  to  the  dimensional  change  along  the  other 
directions,  one  must  be  given  a  distribution  function  of  fiber  orientation,  p(e)  where  e  is 


29 


5.00C-5 


Figure  18.  Comparisort  Between  the  Prediction  of  the  Present  Model  and  the  Experimental 
Data  of  Yoda  et  ai.  for  a  W-Cu  Composite. 


Figure  19.  Comparison  Between  the  Prediction  of  the  Present  Model  and  the  Experimental 
Data  of  Patterson  and  Taya  for  a  SiC/AI  Composite. 
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the  orientation  angle  of  fiber  with  respect  to  the  reference  axis.  Tayaet  al.  (26)  recently 
analyzed  the  dimensional  change  of  a  misoriented  short  fiber  composite,  the  geometry  of 
which  is  shown  in  Fig.  20(a)  where  the  distribution  function  is  assumed  symmetric  with 
respect  to  the  reference  (X3-)  axis,  Fig.  20(b).  The  prediction  of  accumulated  dimensional 
changes  (N  »  1000)  of  a  misoriented  SiC/AI  composite  by  the  model  are  shown  as  a 
function  of  cut-off  angle  p  (see  Fig.  20(b))  in  Fig.  21  where  the  experimental 
measurements  are  also  shown  by  the  dashed  lines.  Despite  the  good  agreement  between 
the  analytical  and  experimental  results  of  the  through-the-thickness  dimensional  change 
(eii),  there  exists  a  poor  agreement  for  the  in-plane  dimensional  changes  (£22  and  £33). 
Namely,  there  is  no  unique  value  of  ^  for  which  the  predictions  of  all  components  of 
dimensional  change  agree  well  with  the  experiment.  This  is  most  likely  because  the 

actual  distribution  of  short  fibers  may  not  be  the  assumed  in-plane  semi-random,  and  the 

symmetrical  distribution  of  short  fibers  about  the  reference  axis  (X3-axis)  may  not  be 
realized  in  the  SiCw/AI  composite.  In  fact,  the  distribution  function  could  be  dependent 
of  both  6  and  the  spherical  polar  angle  (p;  i.e.  a  more  general  distribution  in  three- 

dimensional  space.  It  should  be  noted  in  Fig.  21  that  for  p  >  ff/2,  short  fibers  are 

distributed  in  an  in-plane  completely  random  manner  for  which  both  components  of  the  in¬ 
plane  dimensional  change  coincides,  i.e.  £22  ■=  £33- 

4.2  Dimensional  Change  for  the  Second  Stage 

Regarding  the  second  stage  of  dimensional  change,  an  attempt  is  made  to  compare 
the  predictions  based  on  our  analytical  model  for  progressively  debonded  interfaces  and  to 
limited  experimental  data.^®’^®^ 

Figure  22  shows  a  comparison  between  the  progressively  debonded  interface  model 
and  experimental  results  for  the  W-1%Th02  Fe-25Cr-8AI-0.5Y  composite  of  Fig.  9.^^®^ 
The  experimental  results  were  obtained  for  two  composite  (f  -  .2)  and  two  monolithic  Fe- 
25Cr-8AI-0.5Y  (f  «  0)  specimens.  Here.  Af,^  is  initially  equal  to  zero  but  increases  by  a 
constant  0.0000025  per  thermal  cycle,  while  Sf^  remains  a  constant  0.00002  per  thermal 
cycle.  Therefore,  after  600  thermal  cycles  the  difference  in  effective  fiber  volume 
fraction  between  the  cooling  process  and  the  heating  process  has  increased  to  0.0015,  and 
the  cooling  process  effective  fiber  volume  fraction  has  decreased  from  an  initial  value  of 
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(a)  (b) 

Figure  20.  (a)  Misoriented  Short  Fiber  Composite,  Where  the  Misorientation  Is  In  the  X2- 
X3  Plane,  (b)  Uniform  Distribution  Density  Function  p(e)  Bounded  by  Cut-Off 
Angle  p. 


Figure  21.  Analytical  Results  of  the  Present  Model  for  the  Dimensional  Change  of  a  SiC- 
Al  Composite  After  1000  cycles.  Experimental  Data  of  Patterson  and  Taya  are 
also  shown  for  Comparison. 
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Figure  22.  Comparison  of  analytical  predictions  of  the  extended  model  to  experimental 
results  for  of  the  W-1%Th02  Fe-25Cr-8AI-0.5Y  composite  of  Fig.  6. 


Figure  23.  Comparison  of  analytical  predictions  of  the  extended  model  to  the 
experimental  results  of  Yoda  et  al.  (1978)  for  a  large  number  of  thermal 
cycles. 
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0.200  to  0.188.  Both  of  the  effective  fiber  volume  fraction  quantifications  seem  within 
reason  given  the  damage  porosity  exhibited  in  Fig.  9.  The  analytical  prediction  shown  in 
Fig.  9  assumed  that  damage  development  began  with  the  onset  of  thermal  cycling.  Better 
agreement  between  the  analytical  predictions  and  experimental  results  would  exist, 
however,  if  damage  development  was  assumed  to  begin  after  a  critical  number  of  thermal 
cycles.  Additionally,  the  model  does  not  explain  the  anomalous  compressive  strains 
measured  at  low  numbers  of  thermal  cycles.  The  model  does  however,  correctly  predict 
that  the  strain  accumulated  per  thermal  cycle  becomes  much  larger  as  the  interfacial 
damage  exhibited  in  Fig.  9  becomes  severe. 

Figure  23  shows  a  similar  comparison  between  the  extended  model  predictions  and 
experimental  results  for  the  W-Cu  composite  system.^  The  basic  material  property  and 
thermal  cycle  parameter  inputs  are  the  same  as  that  given  in  Table  4.  For  the  0.10  initial 
fiber  volume  fraction  case,  Afi^was  initially  equal  to  zero  but  increased  at  a  constant 
0.000005  per  thermal  cycle,  while  5fN  remained  a  constant  0.00002  per  thermal  cycle. 
Therefore,  after  1000  thermal  cycles  the  difference  in  effective  fiber  volume  fraction 
between  the  cooling  process  and  the  heating  process  had  increased  to  0.005,  and  the 
cooling  process  effective  fiber  volume  fraction  had  decreased  to  0.08.  Both  effective 
fiber  volume  fraction  quantifications  seem  within  reason  given  the  damage  porosity 
discussed  by  Yoda.^^)  In  the  0.50  initial  fiber  volume  fraction  case,  Af^  was  also  initially 
equal  to  zero  but  increased  at  a  constant  0.00006  per  thermal  cycle,  while  Sf^^  was  a 
constant  0.00015  per  thermal  cycle.  Figure  9  shows  that  the  correlation  between  the 
model  results  and  experimental  measurements  is  very  good  in  the  0.10  initial  fiber 
volume  fraction  case,  however  the  experimental  result  in  the  0.50  initial  fiber  volume 
fraction  case  is  larger  than  that  of  the  model. 

5.  CX>JCLUSIONS 

A  new  thermal  cycling  strain  accumulation  model  was  developed  from  the  previously 
existing  Taya-Mori  modei  (25)  which  gives  quantitative  predictions  for  the  strain 
accumulated  per  thermal  cycle  in  the  case  of  an  aligned  fiber  reinforced  metal  matrix 
composite  with  perfect  interfacial  bonding.  The  present  model  considered  an  initially 
stress  free  elastic/plastic/creep  matrix  reinforced  with  purely  elastic,  aligned  fibers 
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subjected  to  an  idealized  thermal  cycle.  The  present  model  was  found  to  give  results 
similar  to  the  existing  Garmong  (24)  and  Taya-Mori  (25)  models  for  the  specialized  case 
of  continuous  fiber  reinforced  MMCs  thermal  cycled  to  a  maximum  temperature  at  which 
creep  processes  are  rapid.  The  present  model  is,  however,  expected  to  give  predictions 
different  from  the  existing  models  in  the  case  of  a  composite  with  small  fiber  aspect 
ratio  thermal  cycled  to  a  maximom  temperature  where  creep  processes  are  not  rapid.  The 
predictions  of  the  model  were  compared  to  experimental  measurements  of  the  W-Cu 
material  system  and  the  agreement  was  found  to  be  good.  The  extension  of  the  basic 
model  to  the  case  of  misoriented  fiber  reinforcement  was  described,  and  compared  to 
experiment.  Experimental  observations  of  interfacial  damage  development  were 
exhibited,  and  these  observations  were  used  to  formulate  a  progressive  interfacial 
damage  development  process  in  which  interfacial  damage  first  occurs  near  the  fiber  end 
at  a  location  of  maximum  cyclic  shear  strain.  Once  initiated,  the  damage  porosity  quickly 
extends  towards  the  fiber  longitudinal  center.  The  results  of  the  extended  model  show 
that  the  strain  accumulated  per  thermal  cycle  dramatically  increased  from  that  of  the 
perfectly  bonded  case,  and  furthermore,  the  strain  accumulated  per  thermal  cycle 
increases  as  the  interfacial  damage  progresses.  Thus,  the  qualitative  performance  of  the 
extended  model  agreed  with  experimental  strain  measurements  of  thermal  cycled  metal 
matrix  composite  systems. 

Based  on  the  results  of  the  present  study,  the  future  research  that  remains  to  be 
conducted  is  suggested,  as  described  in  Appendix  B. 
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Appendix  A:  Eshelby's  Equivalent  Inclusion  Method 

Here  we  consider  a  composite  which  consists  of  aligned  ellipsoidal  inhomogeneities 
£2  with  a  prescribed  eigenstrain  e^  (called  an  inhomogeneous  inclusion)  embedded  in  an 
infinite  matrix,  Fig.  A(a).  The  stiffness  tensors  of  the  inhomogeneity  (fiber)  and  matrix 
are  denoted  by  CfandCn,.  respectively.  Consider  an  arbitrary  area  in  the  matrix  of  the 
composite  where  the  average  disturbance  of  the  stress,  ‘aOrn  and  strain  e  are  related  to 
each  other  by: 

»  9m  ■  ®  (A1 ) 

Then  introduce  an  additional  inhomogeneous  inlusion  (fiber)  into  .I'e  average  matrix  area. 
The  local  stress  (o)  and  strain  (e)  disturbance  induced  by  the  introduction  of  the 
additional  fiber  is  given  by 

o  =  Cf  •  (e  +  e  -  e^) 

(A2) 

where  e*  is  the  presently  unknown  eigenstrain  fictitiously  introduced  in  the  fiber  domain 
(£2),  which  vanishes  outside  Q,  and  e  is  related  to  e*  by  Eshelby’s  tensor  S 

e-S  e*  (A3) 

Equation  (A2)  indicates  that  the  problem  of  Fig.  A(a)  can  be  converted  to  that  of  Fig.  A(b). 
The  local  stress  o  defined  by  eq.  (A2)  can  represent  the  internal  stress  field  associated 
with  any  fiber  (  such  as  the  existing  one).  The  volume  integral  of  the  internal  stress 
(local  stress)  over  the  entire  composite  domain  (D)  must  vanish: 

J  5  dV  -  0  (A4) 

A  substitution  of  eq.  (A2)  into  (A4)  and  use  of  (A1)  leads  to 

e+f(e-e*)-0  (A5) 
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where  f  is  the  volume  fraction  of  fiber.  By  using  eqs.  (A3)  and  (AS),  eq.  (A2)  can  be 
reduced  to 


[(1  -  f)  (C,-Cm)  •  (S-U  +  C,l  •  e*  =  C,  •  e^ 


(A6) 


Then  e*  is  obtained  as 

e*  =  A-^  Cf  •  (A7) 

where  A“^  is  the  inverse  of  the  6  x  6  matrix  A  defined  by 

A.(1-f)(Cf-Cm)  (S-l^)+C,  (A8) 


and  where  is  identity  matrix. 

It  is  noted  in  the  above  calculations  that  C.  S,  I  matrices  are  6  by  6  whileo  and 
are  6  by  1  column  vectors. 


MT:ldk/dy 
1 1/18/91 
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Appendix  B:  Future  Research 


The  study  that  remains  to  be  done  in  the  near  future  is  stated  in  the  following: 

(1)  Document  the  damage  accumulation  in  W-Th02/FeCrAlY  composite  subjected  to 
combined  creep/thermal  cycling  loadings. 

(2)  Based  on  (1),  construct  a  performance  map  of  the  composite  under 
creep/thermal  cycling  loads.  An  example  of  such  a  map  is  shown  in  Fig.  B. 

(3)  Apply  the  present  methods  of  analysis  and  experimental  techniques  to  Ti-alloy 
matrix  composite  system  where  Tn,ax  and  Tmin  are,  respectively,  the  range  of 
900^0  and  -SOX. 

As  to  (2),  each  iso-damage  line  indicates  a  certain  degree  of  damage,  for  example, 
20%  reduced  stiffness  or  10%  reduced  tensile  strength.  The  intersection  points  with  the 
vertical  axis  represent  the  damage  state  of  MMCs  subjected  to  creep  only,  whereas  those 
with  horizontal  line  are  the  damage  state  induced  by  thermal  cycling  only.  Hence,  the 
damage  state  of  MMCs  induced  by  combined  creep  and  thermal  cycling  loads  is  represented 
by  iso-damage  lines  between  the  intersecting  points. 


Rg.  B.  Performance  map  of  an  MMC  under  creep/thermal  cycling  loadings. 
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In  the  case  of  Ti-alloy  matrix  composite,  its  application  is  most  iikely  to  be  the 
skin  structure  of  high  speed  space  vehicles  such  as  NASP  or  the  High  Speed  Civilian 
Transport  (HSCT).  Thus,  the  appropriate  thermai  cycling  range  for  Ti-alloy  matrix 
composites  would  be  -50  ~  900°C.  With  minor  modification  of  the  present  apparatus,  we 
can  evaluate  the  thermal  cycling  resistance  of  Ti-alloy  matrix  composite  by  subjecting  it 
to  large  numbers  of  thermal  cycling  with  and  without  creep  loading 
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